Antiblockade in Rydberg excitation of an ultracold lattice gas 
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It is shown that the two-step excitation scheme typically used to create an ultracold Rydberg gas 
can be described with an effective two-level rate equation, greatly reducing the complexity of the 
optical Bloch equations. This allows us to solve the many-body problem of interacting cold atoms 
with a Monte Carlo technique. Our results reproduce the Rydberg blockade effect. However, we 
demonstrate that an Autler-Townes double peak structure in the two-step excitation scheme, which 
occurs for moderate pulse lengths as used in the experiment, can give rise to an antiblockade effect. 
It is observable in a lattice gas with regularly spaced atoms. Since the antiblockade effect is robust 
against a large number of lattice defects it should be experimentally realizable with an optical lattice 
created by CO2 lasers. 

PACS numbers: 32.70. Jz,32.80.Rm,34.20.Cf 



The prediction of a blockade effect in the Rydberg ex- 
citation of ultracold gases due to long-range interaction 
has sparked several experiments which have demon- 
strated the blockade effect under different circumstances. 
In Rydberg excitation to n = 30-80 was achieved by a 
single photon of a pulsed UV laser from the 5s Rb ground 
state. In 0, a two-step scheme is used where a strong 
pump laser drives the 5s-5p transition of Rb atoms while 
a tunable second laser excites from 5p to a Rydberg ntt 
state. In all three cases a suppression of the excitation 
has been observed as a function of increasing laser in- 
tensity or/and density of the g cLS 7 1.6., cLS 3j function of 
increasing effective interaction between the atoms which 
shifts the energy levels out of resonance with the laser 
light. 

However, the two different excitation schemes (single 
UV pulse and two-step excitation) have a very different 
nature and may lead to dramatic differences in the block- 
ade behavior. In fact, as we will show below, the two- 
step scheme (see Fig.^l may even induce an antiblockade 
effect due to the form of the single Rydberg atom pop- 
ulation P e (t, A) as a function of the detuning A from 
resonant excitation. This antiblockade effect will only be 
uncovered by a structured ultracold gas, e.g., a lattice gas 
with atoms regularly placed on a 3-dimensional lattice. 
Yet, the condition on the regular arrangement is not too 
strict: even with 20 % lattice defects the antiblockade is 
still visible. 

The key to the antiblockade effect is the excitation dy- 
namics of the final Rydberg state nl in a three-level sys- 
tem with the transition between level \g) (ground state) 
and level \m) (intermediate state) driven by a strong op- 
tical transition with Rabi frequency lu, the Rabi 
frequency which drives the transition between level \m) 
and the Rydberg level |e), see Fig. ^ 

Experimentally, the intermediate level \m) decays with 
a rate r ^> u>, large compared to the upper Rabi fre- 
quency to. Under the conditions Q 3> to and f > w 
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FIG. 1: Sketch of the two-step excitation scheme. 



and a sufficiently long laser pulse the optical Bloch equa- 
tions for transitions in this three-level system reduce 
to a rate equation for a two-level system which contains 
the upper state |e) and an effective lower state \g) [fj. 
The reason for this simplification lies in the strong de- 
cay r which damps out coherences relatively fast and 
ultimately allows for the rate description which directly 
gives the probability P e for an atom to be in the upper 
Rydberg state |e), 



P e (t,A) =Foo(A) 1-exp 



7 (A)t 



Poo (A) 



(1) 



where P^ = P e (t — > 00, A) is the steady-state occupa- 
tion of |e) and 7(A) is the rate of populating |e) for short 
times t. 

Typical Rydberg populations as a function of detuning 
A are shown in the left part of Fig. which also demon- 
strates that, for the given Rabi frequencies and decay 
rate, a pulse length of t > 0.5 /xs is enough to make the 
description with a rate equation applicable. 

For atoms embedded in an environment of other atoms, 
e.g., a gas, the Rydberg level |e) of a ground state atom is 
shifted by S due to the (weak) interaction with Rydberg 
atoms in the neighborhood. 
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FIG. 2: The population P e of the Rydberg state |e) in the three-level system of Fig. according to the rate equation Eq. Q 
(solid) and the optical Bloch equation (dashed) for laser pulse lengths of 0.1/xs (a,e), 0.3/xs (b,f), 0.5/iS (c,g), 1.0/xs (d,h). The 
parameters are tt = 2vr • 4 MHz, lo = 2tt • 0.2 MHz, r = 2tt ■ 6 MHz for the left set (a-d) and Q = 2tt • 22.1 MHz, w = 2vr ■ 0.8 MHz, 
T = 2vr • 6 MHz for the right set (e-h). 



We use the simple picture as formulated in [jj for Rb. 
A pair of Rydberg atoms in state ab at distance R ex- 
periences a shift S(R) of its electronic energy due to an 
induced dipole coupling V(R) — Haa' l^w / R 3 to an ener- 
getically close state a'b'. It is given by the eigenvalues 

5{R) = \{5 ±{5l+W^) (2) 

of the two-state matrix Hamiltonian with elements H\\ = 
0, i?22 = <5o, H12 = H21 = V(R), where So is the asymp- 
totic (R — > 00) difference between the energies of the two 
pairs. 

The relevant neighboring pair to an excited pair ab = 
nsns, corresponding to two atoms in state |e) of Fig. ^ is 
a'b' = (n — ^)P3/2 n P3/2- Eor a specific quantum number 
n we may define [i 2 (n Q ) = Hn s(n ~i)pVn sn p- We have 
taken its value and the value for So (m 2 (ho) = 843800 a.u., 
So = —0.0378 a.u. for n = 48) from [jj and adapted to 
our Rydberg levels by appropriate scaling in n [8j, 

/i 2 (n) = n 2 {n ) ^j^j 

So(n) = S (n ) M?) , (3) 

where n* = n — rj includes the appropriate quantum de- 
fect, for the ns states of Rb r) = 3.13. Furthermore, we 
excite the repulsive branch of Eq. which finally defines 
the shift S resulting from a single excited atom. 

However, we are interested in the many-atom case of 
an ultracold gas, where the laser is tuned to the single 
atom resonance and the detuning A is given by an inter- 
action induced shift 5 of |e). Since the shift is additive it 
is easily generalized to this case where now Si for atom i 
is given by the sum Si = J^j ^(|*t — with the index j 
running over all Rydberg atoms existing at this time and 



5(R) given in Eq. J2J). Number and location of Rydberg 
atoms at a time can be determined by solving the rate 
equation Eq. (|TJ for each atom with a Monte Carlo tech- 
nique. The result is shown in Fig. [3 namely a decreasing 
fraction of Rydberg atoms with increasing level shift S, 
here realized through increasing excitation n. This is the 
typical blockade effect, observed in pi El 
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FIG. 3: The fraction of excited Rydberg atoms f e as a func- 
tion of increasing excitation n and for different laser pulse 
length of I/is (solid), 2^ (dotted), and 5/is (dashed). The 
density of the ultracold gas is p = 8T0 9 cm -3 , and the parame- 
ters are close to those of the experiment Q, & = 2-7T-22.1 MHz, 
lo = 2tt • 0.8 MHz, T = 2tt ■ 6 MHz. 

The parameters of Fig. [3] are close to those of the ex- 
periment jjj. However, for those parameters, the single 
atom excitation probability P e (t,A) differs qualitatively 
from the one shown on the left part of Fig. [2] It has a 
double peak structure due to an Autlcr-Towncs splitting 
induced by the strong driving Q, as can be seen on the 
right part of Fig. with maxima at A = ±A m . Due to 
the wide distribution of mutual atomic distances R in a 
gas the characteristic double peak excitation profile with 



a peak separation of A w f2 does not have an effect on 
f e as shown in Fig. 

To make the effect of the double peak in the excitation 
profile visible, one needs a situation where the distribu- 
tion of mutual atomic distances of the atom is restricted. 
This can be achieved, e.g., with a lattice gas, i.e., with 
atoms regularly spaced on a three dimensional lattice. 
The fraction of excited atoms on a simple cubic lattice 
with a lattice constant a = 5fj,m is shown in Fig. 0] One 
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FIG. 4: The fraction of excited Rydberg atoms f e as a func- 
tion of increasing excitation n for atoms regularly arranged 
in a 3-dimensional simple cubic lattice with lattice constant 
5/xm (p — 8 ■ 10 9 cm -3 ) with parameters as in Fig. [3] (a) per- 
fect filling of lattice sites, (b) 20% lattice defects (i.e., empty 
lattice sites). 

can see the known Rydberg blockade but on top of it an 
antiblockade effect, i.e., an enhanced fraction of excited 
atoms for certain values of n. A closer look at Fig. 0] 
reveals that there are even satellite peaks next to the 
main peaks. This pattern of peaks is easily understood 
if one analyzes the geometry of the atoms placed on a 
cubic lattice (see Fig. 0. The shift in the Rydberg level 
of a ground state atom is dominantly induced by a Ry- 
dberg atom located on the next lattice site at a distance 
of the lattice constant a. Hence we may write for the 
shift 8(a;n). For a certain n of the Rydberg atom this 
shift matches the detuning A m needed for maximum ex- 
citation rate (see right part of Fig. [5J|, which leads to a 
peak in Fig. 01 As experimentally also done, by changing 
the laser intensity we can reach different n keeping the 
Rabi frequency ui constant. Clearly, with different rii the 




FIG. 5: Arrangement of Rydberg and ground state atoms 
leading to antiblockade peaks and satellite peaks. 



(optimum) shift A m can be achieved by a Rydberg atom 
at a different distance <2j, so that in general 



8(a l \n i ) 



(4) 



The obvious candidates in a cubic lattice apart from 
the nearest neighbor at ai = a are the neighbors along 
the diagonal at a 2 = V2a and along the cubic diago- 
nal at 03 = \/3a. If one calculates the corresponding 
quantum numbers Uj from Eq. Q, one would predict 
{ni,n 2 ,n 3 } = {65,78,87}. This differs at most by one 
quantum of excitation from the numerical values of Fig.|3 
which are {65, 77, 86} for the shortest pulse length (1/xs). 
Of course, for the longer pulses the interaction is stronger 
with a considerable amount of excited atoms. This back- 
ground of excited atoms leads to a general shift of the 
peaks towards lower n. 

In addition to the main peaks described, satellite peaks 
can be seen to the left of the main peaks. They corre- 
spond to a situation where two Rydberg atoms on each 
side of a central atom contribute half to the shift of the 
excited level of the central atom, see Fig. 03 Of course, for 
this enhancement less interaction is necessary than in the 
nearest neighbor enhancement. Therefore these satellite 
peaks appear at smaller n, to the left of the main peaks. 
Since double the number of Rydberg atoms is necessary, 
the satellite peaks are roughly only half the height of the 
main peaks in the linear regime of low concentration of 
Rydberg atoms (short laser pulse, solid curve in Fig.|3). 

A lattice with characteristic spacing of a ~ 5/^m can 
be produced by CO2 lasers 0, but probably can only 
be filled leaving a significant amount of lattice defects. 
Therefore, we have simulated a situation with 20% ran- 
dom lattice defects. The defects hardly diminish the 
contrast of the antiblockade effect in the fraction of ex- 
cited atoms (Fig. 0Jd), since the antiblockade relies dom- 
inantly on contributions from pairs of close neighbors: If 
a neighbor is missing due to a lattice defect this does 
neither lead to antiblockade nor to blockade and there- 
fore does not spoil f e . The large tolerance with respect 
to lattice defects makes an experimental verification of 
the antiblockade feasible, depending on the appropriate 
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FIG. 6: Phase diagram for P e (t, A) from Eq. Q. Solid lines 
represent the phase boundaries between the blockade (upper 
area) and the antiblockade (lower area) regime for different 
pulse lengths, while V — 2n ■ 6Mhz. The dashed lines are the 
linear approximations for Q 3> u, see text. The parameter 
sets used in Fig. |5]are marked with black dots. 

combination of Rabi frequencies, pulse length and decay 
rate T. The corresponding "phase diagram" with the 
boundary between the single peak (blockade) and double 
peak (antiblockade) domain is shown in Fig. El The rate 
equation Eq. (JIJ reveals that for large the expo- 

nential is suppressed and the steady state probability P^ 
dominates, giving rise to a single peak structure. On the 
other hand, for small 7i/P>o the Rydberg population is 
governed by 7 which gives rise to a double peak struc- 
ture. The rate 7 as well as P^ take a relatively simple 
form in the limit ui <C T <C Q, namely 



' 2(1 -4A 2 /^ 2 ) 2 00 l + 8A 2 /ft 2 ' w 

The coefficients Eq. J5J reveal a universal condition for 
the transition from the double to the single peak struc- 
ture of P e (f,A), defined by d 2 P e (t, A)/<9A 2 | A =o = 
which can be written as 

,g = 21n(l + . 9o ) (6) 

with go = Ttuj 2 /il 2 . Equation JJJ is easily solved by it- 
eration to give go — 2.513. Hence, for fl ^>> to we expect 
a linear line u> = afl separating the two regimes, where 
a 2 = go/(tr) which is indeed the case (dashed lines in 
Pig. |SJ|. Note also, that Eq. JSJ) clearly demonstrates the 
transient character of the double peak structure which 
vanishes for long laser pulses. Yet, experimentally acces- 
sible parameters, e.g., in realize exactly the transient 
regime and therefore provide the conditions to see the 
antiblockade. 

To summarize, we have derived a rate equation for the 
population of Rydberg states in ultracold gases by a two- 



step excitation process. The rate describes very well the 
structure of the Rydberg excitation in a single atom when 
compared to the exact solution of the Bloch equations 
including a non-trivial transient Autler-Townes splitting 
in the Rydberg population for certain parameters. 

The validity of the rate equation has allowed us to for- 
mulate the many-body excitation dynamics of Rydberg 
states in an ultracold gas without a mean-field approxi- 
mation and for a realistic number of atoms ^2] as a 
random Monte Carlo process. We can reproduce the ob- 
served Rydberg blockade effect observed previously and 
also its effect on the atom counting statistics ^3 as mea- 
sured in 0, but in addition we have identified an an- 
tiblockade effect due to the Autler-Townes splitted Ryd- 
berg population. We predict that this antiblockade effect 
can be seen in an experiment with a gas trapped in an 
optical lattice created by a CO2 laser since the antiblock- 
ade effect is robust even against a large number of lattice 
defects. In the (realistic) limit of a very small upper Rabi 
frequency u> we could show that the formation of the dou- 
ble or single peak structure in the Rydberg population is 
determined by a universal parameter which allows a sim- 
ple navigation in parameter space consisting of the two 
Rabi frequencies, the decay rate of the intermediate level 
and the pulse length, to achieve the desired peak struc- 
ture in the single-atom Rydberg excitation probability. 
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